function [d, t] = getcsg(dat, th, No)
mp = arrayfun(@(x) x.SourceX, th);
d = dat(:, mp==No);
t = th(mp==No);
end